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Abstract 

We establish a scaling limit for autonomous stochastic Newton equations, 
the solutions are often called nonlinear stochastic oscillators, where the non- 
linear drift includes a mean field term of McKean type and the driving noise 
is Gaussian. Uniform convergence in 1? sense is achieved by applying 1?- 
type estimates and the Gronwall Theorem. The approximation is also called 
Smoluchowski-Kramers limit and is a particular averaging technique studied 
by Papanicolaou. It reveals an approximation of diffusions with a mean-field 
contribution in the drift by stochastic nonlinear oscillators with differentiable 
trajectories. 



1 Introduction 

In E. Nelson [18] Brownian motion is constructed as a scaling limit of a family of 
Ornstein-Uhlenbeck position processes which possess differentiable sample paths by 
construction. The so-called "Ornstein-Uhlenbeck theory of Brownian motion" con- 
stitutes a dynamical theory of Brownian motion. In principle, this result goes back 
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to [20] and [14] and also Chandrasekhar [8] studied this kind of limits. The second 
order stochastic Newton equation is represented as a degenerate system of first order 
Ito equations. For more details see Nelson [18] and the references therein. 
Degenerate diffusion processes which satisfy stochastic Newton equations driven by 
Brownian motion have been investigated in many works e.g. [1, 2, 17]. For applications 
in the setting of diffusions on manifold, see [10]. The broader class of Levy processes 
also includes processes with jumps and hence can explain more phenomena. We proved 
a Smoluchowski-Kramers limit for stochastic Newton equations driven by Levy noise 
see [1, 5] in the absence of mean-field terms. In the PhD thesis of Zhang [21] the 
reader finds references concerning generalizations to an infinite dimensional setting. 
Also it should be mentioned that the Smoluchowski-Kramers approximation carries 
over to nonlinear stochastic oscillators driven by Fractional Brownian motion [6, 9]. 
We would like to emphasize that the mathematical techniques and the mode of con- 
vergence in the works mentioned above differ substantially. 

There are studies prior to this work considering the differential approximation in the 
presence of mean field, one of these is [17]. In that paper, a differentiable approxima- 
tion for the stochastic Lienard equation with mean-field of McKean type is achieved. 
In the theory of differential games with a large number of players or in financial math- 
ematics in particular for pricing models for markets with a large number of traders, 
stochastic differential equations with drift of mean- field type appear naturally [11]. 
Here the mean-field term contains no dissipation. In this paper we pursue a scal- 
ing limit as treated in [4, 5, 18] for stochastic Newton equations driven by Brownian 
motion with a nonlinear drift term ((3K), (3 > 0, where K is globally Lipschitz con- 
tinuous, moreover, it depends only on the position process and a mean-field term in 
form of the expectation value of the position process. 

2 Preliminaries 

In this section we recall a representation for Ito processes with linear dissipation and 
explain the scaling limit we are using. 

Proposition 1. Consider the equation 

dY t = -aY t dt + f(t)dt + dB t , 

where a > 0, B t , t > 0, is standard Brownian motion, and f : M — >• M is a measurable 
function such that J T e as \f(s)\ds < oo. Then we know that the solution exists and is 
given by 

Y t = e - at x + [ e- a{t - s) f(s)ds+ [ e- a(t - s) dB s , (1) 
Jo Jo 
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for some initial value xq. For a proof see for example [19]. 

A stochastic Newton equation is an Ito equation which is second order in time, hence 
is given by a system of a first order ordinary differential equation and a stochastic 
differential equation. The dependent variables are traditionally denoted by x t and vt, 
respectively, whence the independent variable t > is interpreted as time. In physical 
models x t describes the position of a particle at time t > 0. It is assumed in this paper 
that the velocity ^| = v exists and satisfies the so called Langevin equation with a 
nonlinear drift depending on the position and its expectation value. 
We introduce a scaling in form of a parameter 7 > for the stochastic Newton 
equation with mean-field having the following form 

dx t = v t dt 

dv t = -jv t dt + K(x t - E [x t \)dt + ^dB t . 

Deviating from Narita, in our model the nonlinear part of the drift includes the mean 
field term E [x t ], moreover, it does not contain the marginal process Vt which is why 
this term does not explicitly contain any scaling parameter. We emphasize that the 
deterministic system corresponding to this system apart from the dissipation 71^ does 
not depend on the scaling parameter. Moreover, for K(x) = — ^ we may define a 
Hamiltonian H(x,v) = \v 2 + V(x), such that ^ = ^ = v and ^ = —(3K(x). 
Averaging properties of such Hamiltonian systems have been studied in [3, 15]. 
For the solution of (2) we scale the time according to tf = -t, and we define the scaled 
process x 7 (t') = x(jt') and f 7 (t') = 71; (7^). Moreover, we introduce a new Brownian 
motion B t > = -^B^t') for a new Brownian motion B. Then (x^ (t') , v 1 (t')) satisfy 
the following stochastic differential equation 

dx], = vj f dt' (2) 
dvj, = -j 2 v] f dt' + ~f 2 K{xJ f - E[x],])dt' + 7 2 ^, 

where in a final step the Langevin equation has been multiplied by 7. For the sake of a 
short notation we replace t' by t, B t , by B t , and set 7 2 = (3. Then the two differential 
equations combine to the initial value problem: 

dx p t = v p t dt (3) 
dvl = -0vldt + pK(x% - i4)dt + 0dB t , (4) 



denotes the expectation, the 



with initial state (xq,Vq) = (xo,Vq), where /if = E 

possibly nonlinear function K satisfies a global Lipschitz condition, i.e. \K(x — y)\< 
k\x — y\ for k > 0, and B t is standard Brownian motion. 
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3 Formulation of the main result 



The system we finally derived will be investigated in this work. When examining the 
equation in the variable Vt, apparently, the drift has the time scale f3t and f3B t = B^ t 
has the faster time scale f3 2 t. Thus, for /3 tending to infinity the time is sent to infinity 
while the Brownian paths performs arbitrarily fast oscillations around the trajectories 
of the deterministic system given by the drift. 

Let us formulate our main result. In a first step we introduce the Ito process in JR 
solving the stochastic differential equation 



dy t = K{y t - E [y t ]) + dB t , 
with y(0) = xq and K as in (3). 

Theorem 1. Let be the solution of (3) and let $ 



(5) 



(xq, vq) be any two- 



dimensional random vector independent of the Brownian motion B t , t > 0, such that 

E [|$| 2 ] < M < oo. (6) 



Then 



E 



sup 

0<t<T 



Vt 



—7-0 as j3 oo, 



for every T < oo, where y(t) is the solution of (5). 

The idea of the proof is to use Gronwall lemma repeatedly to show that \i t and 
E\\ sup 0<t<T Xf] are uniformly bounded on [0, T\. Before we proceed let us make the 
following important remarks. 

Remark 1. As stated in the introduction the function K is supposed to be globally 
Lipschitz continuous throughout the paper. Moreover, assume all conditions stated 
in Theorem 1 hold. Then by Arnold [7], and for more recent results see [12], the 
differential equations (3) and (5) have a pathwise unique solution (x^(t), v^it)^ with 
initial states (xq,v ) satisfying the moment estimates 



E 



sup 

0<t<T 



,0 



< L 



73 > 



for every T < oo with a constant Lp > depending on (3 and T. 

Remark 2. For the same assumptions as in Remark 1 we have for y t defined in (5) 
with y(0) = xq: 



E 



sup 

0<t<T 



Vt\ 



< L, 



for every T < oo and a constant L > 0. 
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Remark 3. The process y t , t > 0, is defined in the same space as the coordinate 
processes xf, t > in (3). In differential geometrical terms this means that the 
driving Brownian motion B t , t > 0, changed from the tangent space M where the 
coordinate process Vt, t > 0, is defined to the manifold which trivially is JR. Moreover, 
the nonlinear vector field K{x t — IE [x t ]) changes from the cotangent space to the 
tangent space. 

Using Proposition 1 the integral equation corresponding to the first equation of (3) 
becomes 



p 



pt pt pS pt pS 

+ v e~P s ds + f3 / e-^"- u) K{x p u - tf)duds + (3 / e~ p{s - u) dB u ds. 
Jo Jo Jo Jo Jo 

(7) 

We change the order in the double integral to obtain 



+ ijj + (3 f e Ps K (x£ - tf) f e-^duds + [ e Ps f e~^ u dudB s , 

Jo Js Jo Js 

where xo = x^(0) and Iq (t) = 1^(0, t) = ^(1 — e _/3t ). Partial integration reveals 

xf = xo + lg{t) + if (t) + [\l - e-^)K{x? - ^)ds + B u (8) 

Jo 

where lf(£) = — e _/3 * J * e l3s dB s . For the sake of a short notation we sometimes drop 
the time parameter, i.e. if := if (•). 



4 Auxiliary estimates 



We evaluate the upper bound for the expectation value of each if (t), % = 0,1, given 
above. For the deterministic integral if the k th absolute moment, k e IV, is trivially 
estimated 



E 



(9) 



Moreover, for if (t) we use Ito isometry to have 











E 


llfWI 2 


< E 













2/3 



(1 



)< 



2/3' 



(10) 
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Using the Lipschitz condition on K we estimate |xf | in (8) by: 



\%t \ — \ x o\ + ~p + \Ii I + 1-^*1 + K / l x = 



By monotonicity of integration there holds \fit\ < E 
and inserting (9) and (10) reveals 



.1 4- 



tf\ds. (11) 
. Taking the expectation 



E 



<E[\x \} + \ + E 



E[\B t \] + 2KjlEl\x?\] ds 



where we used the fact that -2? [|^o| 2 ] < ^ an d -2? < Next we consider the 
moment of the supremum of the square of the solution ar (t), for which we derive a 
bound uniform in the parameter /3 and t G [0,T]. 

Lemma 1. For i/ie same assumption as in Theorem 1, let (x^vfj be the solution 
of (3) with initial state $. Then we have 



sup E 

/3>1 



sup \x (t)\ 

0<t<T 



< H{T) 



for arbitrary fixed T < oo, where H(T) is a positive constant independent of (3 as (3 
tends to infinity. 



Proof of Lemma 1 

Let T < oo be arbitrary but fixed. Consider any t G [0, T]. We return to equation (11), 
square both sides, apply Jensen's inequality before taking the supremum over the 
interval [0,t], and estimate: 

sup |xf | 2 < 5 ( M 2 + tL+ sup |if | 2 + sup \B u \ 2 + 2k 2 T [ \x p s \ 2 + E [\x p s \ 2 ] ds 

0<u<t V P 0<u<t 0<u<t Jo 

(12) 

where we have used that we have positive integrands and the estimate | [it 1 2 < dE \ x t \ 2 

Since B t is a martingale, Doob's inequality reveals E [sup 0<M<t |-B M | 2 ] < 4iE [|£> t | 2 ], 
respectively, by inserting (7) and the recent estimates into 



i/-f/\iO iiO iiO / 

sup | Ji (w) | <2 sup \B U \ +2 sup \B U \ t sup (1 — e 

0<«<t 0<«<i 0<u<t 0<u<t 



-2f3u\ 
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/ E 


sup \x^ 2 


ds^j 


Jo 


0<u<s 





we obtain that for all < t < T: 
E \ sup \x^\ 2 } <$(e \\x \ 2 ] + + 9T + 8T 2 + Ak 2 T [ E 

0<u<t J V P JO 

(13) 

Let us now turn to the core of the proof Lemma 1. By the assumption on the initial 
values x , v we have 

E [\x \ 2 ] < M < oo. 
For /3 > 1 and hence for /3 — > oo on a given interval [0, T] we define the constant 

H (T) := 5M + 5 + 45T + 40T 2 

which combines all additive constants in equation (13) and rewrite 



E 



sup |ar(t)| 

0<«<t 



< 



H (T) + 20k 2 T f E 
Jo 









f E 


1 /3| 2 

sup \x^\ 


ds 


Jo 


0<u<s 





Gronwall's lemma then reveals: 

E 



sup |x^(t)| 

0<u<t 



<H(T), 



(14) 



with H(T) = H (T)e eT2 and 9 := 20k 2 where k is the Lipschitz constant of K. The 
bound holds uniformly for t G [0, T] and for /3 sufficiently large, e.g. f3 > 1, and hence 
for /3 — > oo, which completes the proof. 
Next we give the proof of the main result of this paper. 



5 Proof of Theorem 1 

Under the same assumptions as in Theorem 1. Let (xf,wf) and y t be the solutions 
of (3) and (5) with initial states (xq,Vq) = $ = (x ,v ) and y = x , respectively. 
Combining (5) and (8) we have 



xi - yt 



f K( X P - i4)ds -Ik (y s - fi s ) ds + lg(t) + /f (t) + if (t), 
Jo Jo 



where I?(t),i = 0,1, are as before, l2(t) 



-6t 



f* e l3s K(xP — pL%)ds, and /jl s 



E [y s ]. We exploit that K is Lipschitz continuous with constant k and rearrange the 
arguments in the norm, to get 

|xf - Vt\ < \lg\ + + \ll\ + H \K(x? - /if) - K(y s - ^ s )\ds 

< K Jo \ x s ~ Vs\ds + k f* - (i a \ds + ELo T i (*) 
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By inserting the definition of /if and /i s and by using Jensen's inequality, in particular 



that |/if - n s \ 2 < E 



\ x s ~ Vs\ 



, we estimate further 



\4 - Vt\ 2 < 5 Wt£ |xf - y s \ 2 ds + K 2 tj*E[\xP - y s \ 2 ] ds + J2 \^^\) ' 

Taking first the supremum and then the expectation value yields for < t < T < oo, 
that 



E 



sup \x H u - y u \ 

0<u<t 



< D E 



sup |x£ - yJ 

0<u<s 



(is + 5 ^ 



i=0 



sup 

0<u<t 



with D := 10k 2 T. Direct calculation together with the previous bounds (9) and (10) 
reveals that the last term on the right hand side is uniformly bounded. To this end we 
proceed in the same way as in (12) and (13), in particular we use that K is Lipschitz 
continuous as well as Jensen's inequality to find 



ll(t) < te~ m e 2 ^(K(x p s -^)) 2 ds<K 2 te~ 2 ^ e 2 ^ s -^\ 2 ds 



< 2K 2 te- 2(St / e 2 V\ 2 ds. 



Since the integrand is positive taking the supremum reveals the following estimate: 
sup \ll{u)\ 2 <2^sup 0<lt<t |xf| 2 ^(l- e- 2 ^) < K 2 Tsup 0<u<t \xP\ 2 ± 



0<u<t 



Taking expectation and inserting the bound (14) we obtain: 

„2 









2 


E 


sup 


ll{u) 






0<u<t 







K 



< -TH{T)e eT . 

/ 

For j3 sufficiently large, e.g. j3 > 1, we get the following bound 

2 



(15) 



i=0 





if(«) 


2" 


sup 




_0<«<t 







< - p HT) 



(16) 



with A(T) := 5 (^K 2 TH(T)e eT2 + |j where we combined the bounds (9), (10), and 
(15). By introducing the bound (16) into (5) we arrive at the inequality 



E 



sup K-^l 2 <-A(T) + D [ E 

0<u<t J P Jo 







[ E 


1 P I 2 
sup < - y u \ 


Jo 


_0<u<s 



ds 



which allows to apply Gronwall's lemma a last time, to give 



E sup Lr£ - y 




A(T)e 



,DT 2 



0<u<t 



with constant D as defined in (5) which holds uniformly on the given time interval 
[0, T] and does not depend on (3 for (3 sufficiently large, while -s tends to zero as j3 
tends to infinity. This concludes the proof of Theorem 1 . 
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